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Abstract 

The effects of flavor mixing in electroweak baryogenesis is investigated in a generalized semiclas- 
sical WKB approach. Through calculating the nonadiabatic corrections to the particle currents it 
is shown that extra CP violation sources arise from the off-diagonal part of the equation of motion 
of particles moving inside the bubble wall. This type of mixing-induced source is of the first or- 
der in derivative expansion of the Higgs condensate, but is oscillation suppressed. The numerical 
importance of the mixing-induced source is discussed in the Minimal Supersymmetric Standard 
Model and compared with the source term induced by semiclassical force. It is found that in a large 
parameter space where oscillation suppression is not strong enough, the mixing-induced source can 
dominate over that from the semiclassical force. 
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I. INTRODUCTION 



Electroweak baryogenesis (EWBG) is a promising scenario for explaining baryon number 
asymmetry in the universe, unlike other scenarios valid at grand unification scale it can be 
tested by the upcoming collider experiments. The Standard Model (SM) of particle physics 



due to the too small CP violation 



xperiment 

mm 



and difficulty in generating a strongly first order 



phase transition [sl, ^, Q| can not be a candidate for baryogenesis. However, EWBG is viable 
in many new physics models beyond the SM such as Minimal Supersymmetric Standard 



Model (MSSM) 



11 



12 



H Q, H, y , 0, Q 



3, y, [l^ and the two-Higgs-doublet models 
etc, as in these models there are extra particles contributing to the Higgs potential and also 
new sources of CP violation. 



An efficient way to generate large baryon number asymmetry is to generate it nonlo 
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21| in which the CP violating 



cally, through a charge transportation mechanism 
charges are first generated inside the bubble wall and then get transported into the un- 
broken phases, where the anomalous baryon number violating sphaleron processes are not 
suppressed. In some models such as MSSM and the two-Higgs-doublet model, the bubble 
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231 for the particles giving 



wall of the electroweak phase transition is typically thick 
dominant contributions, and the particle typical Compton wave length A ~ 1/T is much 
shorter than the wall width ~ (10 — 20)/T where T is the critical teinperature. In 
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271, 
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3l| should be 



this regime the validity of semiclassical approach 
justified for single fiavor case, which provides an intuitively simple description by treating 
particle transportation as motion of WKB wave packages. In a slowly moving and CP vi- 
olating Higgs condensate background, the dispersion relation for particles and antiparticles 
are modified differently, contributing to different semiclassical forces, which leads to a net 
excess or deficit of the particle number which can be converted into left-handed fermion 
number asymmetry. The asymmetry of the local fermion density then get transported in 
front of the bubble wall, which bias the baryon number violating processes. 

The problem becomes more involved in the multiple fiavor case, where the CP violating 
mass matrix has nontrivial space-time dependence. The CP violating effects may show up 
not only in the dispersion relations, but also in the mixing of states. In the conventional 



WKB approach 
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3l| the particle propagation is treated as an adiabatic motion: 



the quasiparticles are first rotated into the local mass eigenstates, then these particles are 
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assumed to evolve individually without interference. The semiclassical force terms are iden- 
tified for each individual states as the unique CP violating source for the diffusion equations. 
However, note that in many cases such as in MSSM, the quasiparticles do not decouple even 
in the local mass basis, because the rotation matrices diagonalizing the mass matrix are also 
spatially varying, which leads to interferences among the local mass eigenstates. Namely, 
there always exists nonvanishing off-diagonal terms in the equations of motion. The off- 
diagonal terms in the equation of motion lead to nonconservation of net currents, which 
contributes to extra CP violating source terms. Such a contribution can be numerically 
important, because the semiclassical force term is of the second order in derivative expan- 
sion, and is further suppressed to the third order when entering the semiclassical diffusion 
equations as a source term, while the mixing-induced source presents from the first order. 

The problem of spatially dependent mixing has been discussed in the non-equilibrium 
field theor y appr oach based on the Kadanoff-Baym equations for Wigner transformed Green 



functions 



34 



35| . It has been noticed recently that the evolution of the off-diagonal densities 



from the transverse part of the constraint equation exhibit oscillation behavior, which is 



quite different from the diagonal densities 
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33] . Other methods although using the same 



Schwinger-Keldysh forma! 



condensate insertion 
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ism treat the problem in different manner, such as using the Higgs 
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38| or applying re-summation for Higgs condensate insertion 



but using phenomenological ways to deal with the source terms [39|. The final result still 
differ significantly, especially in the bosonic sector. 

In view of the current theoretical situation, a carefully reanalysis of the semiclassical 
method may still be useful. In this work we generalize the semiclassical method by taking into 
account the non-adiabatic corrections from the spatially varying fiavor mixings. We restrict 
ourselves in the parameter region where the mixing between the local mass eigenstates are 
relatively small, and can be treated as perturbations to the equation of motion for local 
mass eigenstates. Our results show that the correction leads to an extra CP violating source 
which exhibit a typical oscillation behavior in analogy to neutrino oscillations, which is 



similar to the observations made in 



32| . We go a step further to estimate the source term 



and compare it with the semiclassical forces in the context of MSSM. The results show that 
despite the oscillation suppression, the mixing-induced source term has significant numerical 
importance for a large parameter space. 
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II. THE FERMIONIC CASE 



We begin with a brief review of the semiclassical description for fermions in a expanding 
planar bubble wall. For simplicity the fermion is boosted into a wall frame in which the 
momentum parallel to the wall is vanishing. In this case the equation of motion is reduced 
to 1+1 dimensional, and in flavor basis it is given by 

[t{lodt + 73^.) - M\z)Pl - M{z)PR]ij = (1) 

where M{z) is the spatial varying mass matrix from a mass term iPlM{z)iIjr induced by a 
slowly moving Higgs condensate, and Pl,r = (1 =F 75)/2 the chirality projection operator. 
The wave function in the chiral representation of Dirac matrices can be decomposed in terms 
of spin eigenstates 

i^s = e-^"* I ® X. (2) 




with (TsXs = sxs and s = ±1 for two spin states. Ls{Rs) is the left- (right-) handed 
component. Substituting ips into Eq.([T]) and eliminating the right- (left-) handed component 
one obtains the equation for Lg^Rg) in flavor basis 

[uj^ + dl- MM^ + isMd,M-\uj - isd,)] = (3) 
[uj'^ + dl - M^M - tsM^d,M^'\uj + isd,)] Rs = (4) 

One can rotate the flelds into a local mass basis (Lf , i?f ) by defining Lg = UL'^ and Rg = VRg 
in which the mass matrix M{z) is diagonalized WM{z)V = m{z) with m{z) a diagonal 
mass matrix. The phases of U and V are arranged such that the mass eigenvalues rrii are 
real positive. Note that both U and V are z dependent. The equation of motion for local 
mass eigenstates Lf and and Rg becomes 

i{dt - sd, - sU^d,U)Li - mR^ = (5) 
i{dt + sd, + sV^d,V)R'^ - mL'^^ = (6) 

It is evident that the equation for right-handed component in the local mass basis can be 
obtained by replacing U ^ V and s — s. Following the similar step in obtaining Eq.dH) 
one arrives at the decoupled second order equations 

[cu^ + dl -m^ + + isA{uj - isd,)] = (7) 
[a;2 + dl-m^ + 2Ud, - isB{u + isd,)] Rg = (8) 



where 



E = U^d.U , A = U^MdM-^U (9) 
n = V^d.V , B = V^M^dM^-W (10) 

In the above equations the terms of second order in derivative expansion are neglected. The 
off-diagonal elements of A and B are related to S and 11 through 

Aij = — Hij — T^ij , Bij = — Ejj — Uij ( for i ^ j) 

In the remainder of this section we shall suppress the index "d" in the wave function. In two 
flavor mixing case the explicit form of the equation array for left-handed component can be 
rewritten as 

(^"^■^)f'-'Vo (11) 

with 

Du = uj^ + dl -m\ + 2^iid, + isAu{uj - isd,) (12) 
1^12 = 2Si29, + isAuiuj - tsd,) (13) 
D21 = 2S2i9, + isA2iiuj - isd,) (14) 
D22 = uj^ + dl -ml + 2S225, + isA22{oJ - isd,) (15) 

Since the off-diagonal elements -Di2(2i) contain differential operators, one can not obtain 



decoupled equations for Lsi(2) separately. In the conventional approach [25|, the equation 
array is solved approximately by first imposing the derivative operator 1^22(11) on the first 
(second) row of Eq. lfTT]) once again, then picking up the terms up to the first order derivative, 
which leads to equations containing higher order derivatives DaDjjLgj — 0{i 7^ j) and finally 
argue that they reduce to a decoupled form 

DnLsi = 0, (16) 



which is essentially the diagonal part of Eq. pTil . Note that in this method the off diagonal 
terms are counted as higher order terms and thus neglected. An advantage of this method 
is that it naturally allows an adiabatic description of the motion of individual particles. 
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However, It is clear that Eq. flTGll misses the information from off-diagonal term -Di2,2i7 which 
is only valid when |-Di2,2i| ^ \Dii — 022] ■ In large mixing case it is not a good approximation 
to Eq.^. 

To investigate the impact of the off-diagonal terms, in this work we adopt an alternative 
method to solve Eq.(|TTi) approximately. We are interested in the parameter region in which 
the off-diagonal elements is non-negligible but relatively small such that it can be treated as 
perturbations. The validity of the perturbation requires that the mass difference between 
the two mass eigenstates are significantly larger than the mixing term, namely S^m <^ 
|mf — m?j\, where m is the averaged mass. Although this method does not apply to the 
resonance case \mf — m|| ~ 0, it should illustrate the main feature of the mixing caused by 
off-diagonal term Dij. 

Taking the off-diagonal terms D12 and D21 as perturbations. The solutions can be written 
in a generic form 

Ls. = + (17) 
where Lf^ is the lowest order solution satisfying Eq. (fT6ll . i.e. 

DuLf^ = (18) 

and are the corrections due to the off-diagonal terms. The lowest order solution L^°^ is 
obtained by the usual WKB wave ansatz 

where pd is the canonical momentum and the function Wi provides the correct normalization 
for L^^. The real and imaginary part of Eq.ffTSl) lead to two separated equations 

22 2 



- Pli - = Im [2pciT.ii + s{uj + spci)Aii] (19) 



w: 



p^ + 2pe^— = Re[2p^^J:ii + s{u + spci)Aii] (20) 



where the notation prime stands for the spatial derivative dz- The first equation gives the 
dispersion relation 

S ( Cc^ j f^T^ 1 

Pci -Po 7: —IraAii - ImEii + a' (21) 

2po [—1 

where a is an arbitrary phase factor from gauge invariance [27||. From the dispersion relation 



one can deduce the group velocity Vg = {duj/dpc)z and semiclassical force Fi = uj{dvg/dt) 



Poi mjlmAi, 
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Note that only the second term in the force term is CP violating, and is proportional to the 
spin of the local mass states. Together with the lowest order solution to Wj from the second 
equation, the wave function at the lowest order is given by 

Lf = ^ e^rP..^^\ and i?f = , e-r^"'^-' (24) 

with = cj^ — . 

Substituting the off-diagonal terms in to the equation, the first order perturbation takes 
the following form 

L,.^lL°' + L« = l2'+.4? (26) 
fl.. ^fl(»)+fl2' = fl2'+5.iL°' (26) 

which are mixtures of the two unperturbed states. The mixing parameters for particle i to 
the first order of derivative are given by 

^ ■ '^^ijPoj + Aji-sio + Poj) and 5^ = i ^^'i^'^i + Bij{-su+pQj) 
* ml — rrtj ' ' mf — m| 

The mixing parameters for particle j can be simply obtained by replacing i j from 
the above expressions. It is clear that the expansion is valid for dz'mi/{Am^) <^ 1 as 
expected. The momentum dependencies comes from the differentiation operators in off- 
diagonal element D12. 

During the particle propagation inside the bubble wall, the Higgs condensate background 
causes mixing between the two local mass states, which keeps the total particle current 
conserved but leads to nonconservation of the currents for individual particles. The nonvan- 
ishing divergence of the current caused by the off-diagonal part of spatially varying rotation 
matrix is proportional to the mixing parameters and 5i, which are different for particle and 
antiparticle when the mixing contains CP violation. This provides an extra CP violating 
source other than the semiclassical force. Note that the mixing term modifies the disper- 
sion relation as well. As it contains derivative terms, it implies the breaking down of the 
semiclassical picture. However, this is a next to leading order correction to the dispersion 
relation and is expected to be small. 
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As an illustration, we calculate the left-handed current divergence. The temporal and 
spatial part of the left-handed current are j° = L*Ls and = —sL*Ls respectively. The 
mixing parameters for the antiparticle can be obtained by replacing J] T,* , A —>■ A* ,11 
n*, and B B* . The mixing-induced left-handed current, after subtracting the antiparticle 
component is given by 

mi-rrij ^ pQipQj[uj + spQi)[uj + spQj) J 
The CP violating force term, as it is proportional to spin s, only contribute to the 
spin-weighted density. To facilitate the comparison with the force term, we calculate the 
spin- weighted mixing-induced source term. Indeed, it has been pointed out that only spin- 



weighted current is sourced by the moving wall 
source term for the left-handed current is 



32| . In the wall frame, the corresponding 



Sl^ = E '-d,3Z (29) 



2mimj 

2 2 
mf — m- 



rrii 



ImEgLipoijPoj) lmIlijgR{poi,poj 



(30) 



with 



■{Poi - Poj) sin J {pcj - Pci)dz' (31) 

9l R[Poi,Poj) = —j^^^=^^^=^^=^^= -j^^^^^^=^^^^^^=- (32) 

sJpQiP^jiy) + PQi){uJ + poi) VPoiPojii^ - Poi)i^^ - Poj) 

which is a momentum odd function. 

The numerical importance of the source term comes from the fact that it is of the first 
order in derivative expansion, larger than that from semiclassical force which is of the second 
order. Furthermore, The final form of the semiclassical force term entering the semiclassical 
diffusion equation is of the third order in derivative as it is weighted by group velocity. While 
in the diffusion equation, the mixing-induced source term is not suppressed. Therefore in a 
naive counting, one expects a much larger contribution from mixing induced source term. 

What might significantly suppress the mixing-induced source is the oscillation. A crucial 
feature is that the oscillation frequency depends on the momentum dijference rather than the 
momentum themselves, which makes the oscillation suppression less effective for relatively 
small momentum differences. From Eq. (l3T]l . the oscillation term can be approximated by 
sin(po2 —Poi)dz at the lowest order. In a highly relativistic limit the oscillation wave length 
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^osc 

is given by 

1 2 2 

1 mf — mi 

~ ^ £ 

Xosc 47ra; 

which is well-known for neutrino oscillations. For an illustration, taking cu ~ 300 GeV, 
mi ~ 200GeV and m2 — 100 GeV, we find Xosc is around C(0.1GeV~^) which is much larger 
than the typical Compton wave length 1/T. Thus the oscillation is unlikely to completely 
erase the mixing-induced source. 

Finally, The source is proportional to the off-diagonal elements of Ej^ and Uij, which is 
different from the CP violating source term which is always proportional to the diagonal 
elements An. This leads to a different parameter dependence in the final source term. 



III. THE BOSONIC CASE 

The mixing in bosons can be calculated in a similar way. The Klein-Gordon equation in 
flavor basis is 

8^^81,(1) + M\z)(l) = (33) 

where M'^(z) is the mass square matrix which is diagonalized by a unitary transformation 
U^M^{z)U — m'^{z). Rotating (f) into the local mass basis with = U4>'^, the equation of 
motion for 0^ in the wall frame has the form 

a^c/.'^ - 8l<P'^ + 2U^d'Ud,(f)'^ + m^(t)'^ = (34) 

We shall suppressed the index "d" in the remainder of this section. In the two flavor mixing 
case the above equation can be rewritten in a matrix form 




(35) 



The operators Dij are given by 





= 8'^ + ml 


+ 2Ena, 














D22 


= 8"^ + ml 


+ 2E22a, 



(36) 
(37) 
(38) 
(39) 
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with S = WdzU. Again we solve the equation by taking off-diagonal elements -Di2(2i) as 
perturbations. The solution takes the form 



+ €' (40) 



with (pf'^ satisfying the lowest order equations 



A.0f ^ = (41) 

Making use of the WKB ansatz 

0f) = u^.e-*^*e^^'P^^('')"'^' (42) 



and substituting it into the Eq. (llT]) . the real and imaginary part gives two independent 
equations 

-uj^ + pI + m\ - ip'^^ + 2ImSii ■ Pd = (43) 
w' 

+ ^—Pci + 2ReSii -pci = (44) 

The first equation indicates the modification of dispersion relation. An important difference 
from the fermionic case is that the dispersion relation for scalar fields will not generate the 
semiclassical force at the first order in derivative expansion because the contribution from 
ImS cancels in the expression for dvg/dt. The second equation determines the form of w. 
To the lowest order Wi ~ Taking the mixing term Dij as perturbations, the solution 

takes the form 

0,^0f)+e,0f , (45) 
the mixing coefficient to the first order in derivative is 



2 2 

mf — m^- 



(46) 



The mixing parameter for particle j can be obtained by replacing i j in the above 
expression. The spatial component for the current is jf = i{(j)*d^(pi — (?^0*0i). Substitute 
the expression of (pi the mixing-induced current is given by 

jT = ^{e4'^*^<f>f - e*d<f>f><f>f^) + z{e* <l>f> d<j>f^ - 6*90^ )>f ) (47) 

The divergence of the current is 

= 9,jr = -2lm[e[{d<f>f><pf + >f )] (48) 
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In obtaining the above expression, we have used the lowest order equation of motion for 
(l)f\ Note that the terms proportional to all canceled out. Only the e[ terms remain. 
Therefore the source term is of second order in derivative, which is of the next to leading 
order compared with the fermionic case. This makes the contribution from boson mixing to 
be expected less important. 



IV. SOURCE TERMS IN DIFFUSION EQUATION 

At the semiclassical level the divergence of the currents derived in the previous section 
can be directly added to the Boltzmann equation as semiclassical sources. In the wall frame 
the Boltzmann equation has the form 

where fi = fi(t,x,p) is the phase space distribution function, C is the collision term and the 
Pz the kinetic momentum defined by Pz = ujVg. In a fiuid ansatz, the distribution function 
fi in the wall frame is 

f^= ePbA^+.J.)-M^)]±l +^f^('^^P) (50) 
where the spatially varying chemical potential ^i{z) describe the departure from chemical, 
and the perturbation 6fi{z,p) describe the response to the semiclassical force, which by 
definition has no contribution to the particle density, i.e. fcPp6fi = 0. The factor v^Pz 
comes from the Lorentz boost from the plasma frame to the wall frame. Substituting the 
distribution function into the Boltzmann equation 

^ (v^F.z - ^'^) + ^5/; = C[f„ /, ■■■] + S, (51) 

and boosting back into the plasma frame by a Galilean transformation Vg ^ Vg + which 
is valid for small <^ 1, we have 

+ v^) ( + 5/;) + v.^6Fz. = Cf%, 6f.] + Sf (52) 

with /oi the up-perturbed Fermi-Dirac distribution /o = (e'^'^ + l)"^ in the plasma frame. The 
quantity /ij and 5Fzi contain only CP violating part of chemical potential and semiclassical 
force. From Eq. (l23l) . the CP violating force term is 

{mflmAii)' 

oFzi = -s — (53) 
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The source term Sf in the plasma frame is obtained from Si by a boost Pz Pz + v^u, 
namely Sf = Si{pz + VwUj). Note that the p^— even part in Sf only comes from f^a; term 
which is nonvanishing after momentum integration. Thus the source is proportional to the 
wall velocity Vw and is vanishing in the limit of ft„ — > 0, which is a physical requirement. 
Integrating over dP'p weighted by 1 and Pz/^jJ respectively, we get two transportation equa- 
tions 



T 



Pz 



UJ 



LO I I \ UJ 



where the definitions for the integration are 



27 



3l| 



cn + (sf 



PlCf) + /tlsf 

UJ I \uj 



Pz 



(X) = ki \ for X 

\ / ^ r n df(H ^ 



4,-5i^..), and {X)^k 
uj^ UJ T j <Pp^ 



(54) 
(55) 



U^pX , ^ , , 

^ ( X = others) 



Jd'p'-t 

The collision terms can be expressed in terms of chemical potentials and inelastic interaction 
rates 



With the factor ki = J d?p^^/ J d^p ^^°^^^' ^' which is 1(2) for massless fermion(boson). 



a 



<pi 



(fc) 



Pz ^pl 



UJ 



h{{pz/uj)5U)T\ 



where Ff^ is rate for the inelastic channel {k) and ^j'^'* are the signed sum over the relevant 
chemical potential. The quantity V\ is the total interaction rate. Eliminate the quantity 
i^^fi) differentiating the second equation in Eq. (l55ll once again, one arrives at the usual 
form of the diffusion equation 



Sf = ^^'"wDi /El^p \ Q^^fi g = L /gpi 



with ti = k,Ti 



where ^ = ^^i/T the rescaled chemical potential and Di = {{pz/uj^) /F* the diffusion con- 
stant. In deriving the diffusion we have neglected the {(p.HSf)' term whieh is subleading 
compared with (sf). 
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V. MIXING INDUCED-SOURCE IN MSSM 



In MSSM, the chargino transportation provides a dominant CP violating source to elec- 
troweak baryogenesis. The asymmetry in chargino number is converted into the asymmetry 
in left-handed top quarks through Yukawa interactions. The two Higgsino SU{2) doublets 
are hi = {h\j^,h^^)^ and h2 = {h^i^h^i)^ respectively. Together with the two charged 
gauginos and W£ , the charginos are combined into two left- and right-handed four 
component spinors as ■0^, = i hthi^ ^^id %1)r — {{Wj^Y, (^ll)'^)'^- The chargino mass 
term has the form '4)rM{z)%1)l in the wino-higgsino space with mass matrix 

/ 



M{z) = 



(56) 



M2 gH2[z) 
\gHi{z) n 

where M2 and // are soft supersymmetry breaking parameters containing CP phases and 
Hi{z) and H2{z) the Higgs vacuum expectation values (VEVs). The mass matrix is di- 
agonalized by a bi-unitary transformation V^MU. The explicit form of the two rotation 
matrices are 



^^^^=\ (57) 




i(A + A) -a 



a 



|(A + A) 



2 



^ - I + I (58) 



VA(A + A) I a* i(A + A) 



with 



a — M2U1 + ii*U2 , a — M2U2 + ^iUl (59) 
A = iMsl^ - + ul-ul , A=\M2f - -ul + ul (60) 
A = ^A2 + 4|a|' , Ui = gH,{z) (61) 

The two mass eigenvalues are given are interested 

in the states which is Higgsino-like since only Higgsino can be efficiently transported. The 
Higgsino-like state is identified as particle 2(1) for parameter region M2 > |/x|(M2 < |/i|). 
From the above expression we obtain the off-diagonal elements of matrices E and U which 
is relevant to the mixing induced source. 

S12 = -X(aTa) + " + ^'^"^ ^^^^ 
ni2 = -X(aTA) + " + ^^^^ 
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the other off-diagonal elements S2i(n2i) can be easily obtained as S(n) is anti-Hermitian. 
The quantity An relevant to the CP violating force term is given by 



. M2/i, 
^11 = -^[UlU2> 



For simplicity, we only consider the chargino contribution as a dominant source, as the 
transportation for neutrilinos is suppressed due to much weaker couplings to fermions. The 
stop contributions are also suppressed as they are bosons. Furthermore, the mass difference 
for the two stops has to be very large. The left-handed stop must be above 1 TeV to 



evade the LEP experiment constraints 



while the right-handed stop must be light a 



40 



4l|. The large 



round 100 GeV in order to generate strongly first order phase transition 
mass difference between and leads to a very fast oscillation which strongly suppressed 
the mixing-induced source in stop sector. As the semiclassical force term does not present 
in the first order in derivative, we thus neglect the stop contributions. 
For numerical illustrations we take the following simple wall profile 



uiz)=g^/HKz) + Hliz)=g^-- 



1 ■ z 
tanh . 



with Hi(z) = u(z) sin P and H2{z) = u{z) cos P and g the weak gauge coupling. Vc is the 
Higgs VEV at the critical temperature T which is normalized at vt=o — 246 GeV. For 
numerical calculation we take the following reference values 

Vc = 120GeV, T = 90GeV, tan /? = 3, f ^ = 0.03 

We consider different sets of wall-width and MSSM parameters which is sensitive to the 
mixing-induced source. In FigH], we give the source term rescaled by 1/T from both semi- 
classical force and mixing-induced source term for wall width = 10/T and 15/T respec- 
tively. The MSSM softbreaking parameters are fixed at M2 = 200GeV and = lOOGeV 
corresponding to a chargino mass differences — = 111 GeV in the broken phase. The 
CP phase 0^ is set to a typical of 0^ = 0.02 which is a typically allowed value after the 



36|, Q, Q- Note 



constraints from the electron and Hg electric dipole moments (EDMs) 
that the EDM constraints on CP phases depends on MSSM parameters such as sfermion 
masses and CP-odd Higgs masses. For heavy left-handed stop mass above a few TeV as 
also favored by the LEP limit, the allowed CP phase can be much larger. Furthermore the 
two-loop corrections to EDM are controlled by CP-odd Higgs mass, which can be weakened 
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z z 



Figure 1: (a) Left, the semiclassical force induced source term as a function of z (GeV~ ). (b) 
Right, the mixing-induced source term as a function of z. The two curves in each plot corresponds 
to Lw = 10/T (sohd) and 15/T (dashed) respectively. The MSSM parameters are fixed at M2 = 150 
GeV and /i = 100 GeV, with c/)^ = 0.02. 

by heavy Higgs mass, and has no significant effects on the source term. Therefore, although 
small 0^ is used in the calculation, the possibility of large CP phase close to the maximum 
may still be allowed in this case. The curves given in FigH] show that both of source terms 
have nontrivial spatial dependencies, but their origin are quite different. The variation of 
the force term comes from the third derivative of the kink-type Higgs condensate, which 
has one minimal and two maximums. While the mixing-induced source term varies due to 
both the wall profile variation and the oscillation. The oscillation leads to multiple local 
minimum appearing in the curve and is suppressed at large distance by the wall profile. For 
M2 around 200 GeV, the mixing induced source term peaks at z ~ 0.03 with an amplitude 
Sm/T ~ 1.6 X 10^^, much larger than that from semiclassical force which peaks at 2; ~ 0.08 
with Sf/T = 0.6 X 10-^ 

In Figill we give the results for a larger M2 = 250 GeV, and still fix = lOOGeV, 
corresponding to a larger chargino mass differences mi — = 157 GeV in the broken 
phase. One sees that for a larger chargino mass difference both the source terms becomes 
smaller as they are 1/A suppressed. The oscillation in the mixing-induced source term 
becomes obvious and the wave length is shorter, thus the oscillation suppression is stronger. 
The mixing induced source term peaks at 2; ^ 0.03 with an amplitude Sm/T ~ 8 x 10^^. 
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Figure 2: (a)Left, the semiclassical force induced source term as a function of z(GeV~^). (b) Right, 
the mixing-induced source term as a function of z. The two curves in each plot corresponds to 
Lw = 10/T (sohd) and 15/T (dashed) respectively. The MSSM parameters are fixed at M2 = 200 
GeV and /i = 100 GeV, with c/)^ = 0.02. 

M2iGeV) 200 250 300 350 450 500 
5m(Lu, = lO/T) 27.5 8.29 3.37 1.50 0.243 0.0364 
SpiL^ = 10/T) 1.89 1.23 0.822 0.575 0.317 0.246 
= 15/r) 10.9 5.48 1.45 0.604 0.121 0.0177 
Sf{L^ = 15/r) 0.561 0.365 0.243 0.170 0.094 0.073 

Table I: Averaged source term (in unit of 10~^) for mixing induced source and the force term. For 
wall width = 10/T and 15/r respectively. 

Although significantly reduced, it still much larger than that from semiclassical force v^^hich 
peaks at z ^ 0.08 with Sp/T ^ 4 x 10^^ 

To estimate the oscillation suppression effects it is useful to define an averaged source 
over the wall width 

Sakf) = 7fi- I SM{F){z)dz (64) 

J^w Jo 

We calculate the averaged source numerically for different M2 = 200 ~ 500 GeV and list the 
results in TabUl One sees that Sm dominates over Sp in the range 200GeV < M2 < 350GeV. 
With the value of M2 increasing, the averaged source term Sm drops rapidly. For = 10/T, 
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at M2 — 350 GeV, the mixing induced source is only 5.4% of that at M2 = 200 GeV, the 
suppression is due to the increased oscillation frequency. The suppression in force term Sp 
is mainly from the factor 1/A, which decreases much slower. At M2 = 350 GeV, it is still 
about 30% as large as that at M2 = 200 GeV. At M2 — 200 GeV, their relative size is 
Sm/Sf — 14.6. For large M2 — 350 GeV, although they are close in size, the mixing source 
still dominates with Sm/Sf — 2.63. This dominance has a mild dependence on the wall 
width. For — 15/7, the relative size between the two kind of sources remains roughly 
the same, although both of the source term becomes smaller. The mixing-induced source 
(semiclassical force) term is ~ 40%(~ 30%) of that at = 10/T. When the value of M2 is 
around 450 GeV, the two type of source term becomes comparable in size. For a very large 
M2 — 500 GeV, the semiclassical force term becomes dominate, and the mixing-induced 
source term is about an order of magnitude smaller. 

In conclusion, we have studied the effects of flavor mixing in a generalized WKB ap- 
proach in which the off-diagonal terms in the equation of motion are taken into account as 
perturbations. With the presence of a slowly moving CP violating bubble wall, an extra 
mixing-induced CP violating source appears which exhibit an oscillation behavior in analogy 
to the neutrino mixings. In two flavor mixing case, the oscillation frequency is proportional 
to the difference of mass square. The size of the mixing-induced source could be larger than 
that from the conventional semiclassical force as it is at the flrst order in derivative expan- 
sion. We have made a numerical study of the oscillation suppression effects for chargino 
case in MSSM. For a small = 100 GeV, in a large range 200 < M2 < 350 GeV, the mix- 
ing -induced source term dominates over the semiclassical force term. The method is valid 
for small mixing case, a more general method is needed to deal with the maximum mixing 
case. However, even in the small mixing case the mixing-induced source already indicate 
that a signiflcant enhancement of the flnal baryon number asymmetry is possible, which will 
relax the tension between the observed baryon number asymmetry in the universe and the 
constraints from electron and atom EDM in the MSSM. 
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